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I— I Abstract 

Proton Compton scattering is analyzed with the chiral Lagrangian. Partial- wave am- 
plitudes are obtained by an analytic extrapolation of subthreshold reaction ampli- 
(— I tudes computed in chiral perturbation theory, where the constraints set by electro- 

'— ' magnetic-gauge invariance, causality and unitarity are used to stabilize the extrap- 

olation. We present and discuss predictions for various spin observables and polar- 
izabilties of the proton. While for the transition polarizabilities Jeim2,Jmie2 we 

^ recover the results of strict chiral perturbation theory, for the diagonal "Jeiei ■, Jmimi 

elements we find significant effects from rescattering. 

m 

O 

^ 1 Introduction 

> 

^ Nucleoli Compton scattering serves as an important test for Chiral Perturba- 

tion Theory (^PT). One of the reasons is that for such a process the leading 
order loop contributions do not depend on unknown parameters, when cal- 
culated in the heavy baryon formalism [Ij. The same is true for some of the 
higher order terms [2|3] . The application of strict xPT is however limited to 
the near threshold region. A method to extrapolate xPT results beyond the 
threshold region using analyticity and unitarity constraints was proposed re- 
cently in IHEin]- The results for the differential cross section and beam asym- 
metry for the proton Compton scattering from threshold up to ^/s ^ 1300 
MeV are in agreement with empirical data even though all the parameters 
were determined from vrA^ elastic scattering and pion photoproduction. The 
spin-independent polarizabilities of the proton in this approach do not differ 
from xPT prediction by construction [1] and are equal to ap = 13.0 ■ 10""^ fm^. 
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/3p = 1.3 ■ 10~^fm^ in the heavy baryon formahsm at the order Q^. The spin- 
dependent polarizabihties and further polarization observables such as double 
polarization observables Tj2x and Tj2z were not discussed in j^j. Double polar- 
ization observables are of particular interest since they are best suited to access 
the spin-dependent polarizabihties. Measurements are under way at MAMI [7] 
in the energy range above pion-production threshold, and experiments at lower 
energies are also planned at the HI7S facility |H]. 

It is the purpose of this work to derive such quantities in our unitary and 
effective field theory and compare with previous phenomenological approaches 
based on fixed-t dispersion relations (DRs) PITU] . 



2 Chiral symmetry, analyticity and unitarity 



Our approach is based on the chiral Lagrangian involving pion, nucleon and 
photon fields [TT|T2] . The terms relevant at the order for pion elastic scat- 
tering, pion photoproduction and nucleon Compton scattering are taken into 
account H 

Ant = - ^ ^^7^^ (r • (tt X (9^7?))) iV + iV75 7^ (f ■ {d,^)) N 
- e { (vf X (9^7f )) ^ + iV 7^ iV - iV 75 7^ (r X Tf) ^ iv} 



mlN{7i-n)N- , N {d^ tt) ■ (9,7?) (9^9"iV) + h.c. 



di + d2 



N (r ■ ((9^7f) X (d^d^n))) {d'N) + h.c. 



+ N (f ■ ((9^7f) X (d^dxn))) {d''d>'d'N) + h.c. 

2 J 

- 2 1 iV (f ■ (tt X (9^7f)) ) (d'^N) + h.c. 



P rriN 

^-^f" N [ds {d^ 773) + (r ■ (9,7?))) {d^ N) F^, + h. 



I e 



frriN 

+ ^ N a^'' UdJ) ■ (9^9,7f )) (9^iV) + h.c. 



1 Note a typo in Eq. (1) of [IJ. 
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'-N^,Y(rx {dx if)) (d^d^N) F^, + h.c. . 
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6(^20 
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A strict chiral expansion of the various scattering amplitudes to the order 
includes tree-level graphs, loop diagrams, and counter terms. The latter 
depend on a few low energy constants "Cj' and "dj" , which are to be adjusted to 
the empirical data set. Our goal is to improve the convergence of strict xPT by 
an analytic continuation of the subthreshold amplitudes. As compared to strict 
xPT the applicability domain is extended towards higher energies beyond 
the threshold region [4J. The extrapolation of the amplitudes is performed 
insisting on basic principles of analyticity and unitarity. We seek to separate 
left- and right-hand singularities of suitably chosen partial-wave amplitudes 
with angular momentum J, parity P and channel quantum numbers a, h (in 
our case vriV and 7iV channels). The separation is achieved by means of the 
partial-wave dispersion relation 



where the generalized potential, f/^^ (a/s), is the part of the amplitude that 
contains left-hand cuts only. The phase-space matrix Pctf (w) reflects our par- 
ticular convention for the partial-wave amplitudes, that are free of kinematic 
constraints. The dispersion relation ([2]) relies on the fact that the discontinuity 
along the right-hand cut is given by the unitarity condition 



It is emphasized that the separation ([2]) is gauge invariant. This follows since 
both contributions in (|2]) are strictly on-shell and characterized by distinct 
analytic properties. The relevance of the subtraction scale \im in ([2]) will be 
discussed below. 

In the present approach we take into account only the vrA^ intermediate states 
in ([2]) (channels c, d in the sum). The higher mass intermediate states are 
then effectively included in the potential U^ah'^^\f^)- neglect the 7iV 




(2) 




(3) 
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intermediate state, since the effects thereof are suppressed by the square of 
the electromagnetic charge. 

For a given approximative generahzed potential the nonlinear integral equa- 
tion (|2| yields scattering amplitudes, that obey the constraints set by the 
unitarity condition. On the other hand the crossing symmetry constraint is 
not automatically satisfied. If crossing symmetry would be implemented ex- 
actly, a crossed reaction can be computed from the direct reaction amplitude. 
Necessarily, the crossed reaction would be determined by the direct reaction 
amplitude evaluated at energies below the s-channel threshold. Now, imag- 
ine that in a first step we have some means to approximate the generalized 
potential systematically and very accurately at energies above the s-channel 
threshold only. We may use this potential and determine via ^ the scattering 
amplitude in the physical region. In this case we are not able to compute the 
crossed reaction, since our approximation for the generalized potential is valid 
above threshold only, by assumption. The crossed reaction may be computed 
by setting up the analogous nonlinear integral equation (|2| for the crossed re- 
action. It involves potentials distinct from the potentials of the direct reaction. 
But where is the crossing symmetry constraint? Our argument illustrates that 
crossing symmetry is a constraint that affects dominantly amplitudes at sub- 
threshold energies. In a second step we now imagine that we have some means 
to approximate the generalized potential systematically at energies a/s > /xm, 
i.e. we include a small subthreshold region. In this case crossing symmetry 
does provide a constraint. The crossed reaction amplitude can be computed 
from the direct amplitudes, however, only in a specific subthreshold region. 
The coincidence of the crossing transformed amplitudes and the amplitudes 
from the crossed reaction in that specific subthreshold region defines the de- 
sired constraint. Note that a nonempty coincidence region requires fiM to be 
distinct from both s- and w-channel unitarity branch points. Our particular 
choice fj,M = misf implies a significant size of the coincidence region. 

The crucial assumption of our approach is the perturbative nature of the 
reaction amplitudes at subthreshold energies, where the crossing symmetry 
constraints are relevant. To this extent we satisfy crossing symmetry in a per- 
turbative manner. In order to arrive at even improved coincidence properties 
we insist on a subtraction at ^/s = /^a/ in ([2]). Since the generahzed potential 
of the direct and crossed reactions are evaluated in a perturbative manner, 
the crossing symmetry properties of perturbation theory are transported into 
our approach. 

We recall the specifics of how to approximate the generalized potential |1] . In 
an initial step the various reaction amplitudes are evaluated in strict xPT to 
a given order, in our case Q^. From this result we may extract the general- 
ized potential accurate to the same chiral order. In a second step we apply 
conformal mapping techniques to extrapolate the potential into a larger en- 
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ergy domain Q. This is done by means of a transformation ^(a/s) that maps 
Q onto the unit circle. The domain Q includes energies significantly above 
the s-channel threshold, but excludes energies below the w-channel unitarity 
threshold necessarily. While it is possible to extrapolate the potential towards 
higher energies, it is not straightforward to extrapolate it to energies below the 
M-channel unitarity threshold. This is a consequence of the analytic structure 
of the potential. Since the potential is characterized by left-hand cut struc- 
tures it is an analytic function for energies larger than the s-channel unitarity 
threshold. The presence of all left-hand cut structures imply that the potential 
is not analytic for energies sufficiently below the s-channel unitarity threshold. 

All together the generalized potential can be represented as a series in ^ con- 
vergent in the domain Q with 



U {y/s ) = ?7inside(yS ) + f/outsidc(v^ ) , (4) 
^^outsidelVs ) = J2Uk k{\/s )] , Uk = ^ ^""tl'^'jll ^^^^ 

where we allow for an explicit treatment of cut structures that are inside the 
domain fl. These are contained in f/insidc(A/^) correspond typically to 
contributions from tree-level diagrams. The coefficients Uk are uniquely deter- 
mined from the Taylor expansion of the generalized potential in the vicinity 
of the s-channel threshold. To the order we are working at, we truncate the 
series in Q at the zeroth or the first term depending on the channel and the 
partial wave. 

We close our brief primer on the essential ingredients of the novel scheme intro- 
duced in [3j with a discussion of so called CDD poles (TSUI]. In general there 
is no guarantee that a solution of the non-linear integral equation ^ exists 
and/or is unique. Different solutions may be characterized by the presence of 
distinct CDD poles in the N/D ansatz that solves ([2]). We need a condition 
which singles out the unique and physical solution we are interested in. We 
require that the non-linear equation can be solved in perturbation theory at 
least close to the matching scale fiM- To be specific, the reaction amplitude 
must be accessible via an iterative solution of (|2| in a small domain around 
fiM- Two remarks are here in order. First, such a condition implies that the 
reaction amplitudes are consistent with the crossing symmetry constraint, as 
discussed in detail above. Second, this condition excludes the presence of a 
CDD pole. In a case where this iterative requirement cannot be met we admit 
the presence of CDD poles [T3fT^ . As shown in |1] the perturbative nature of 
the reaction amplitudes can be restored by a suitable correlation of the CDD 
pole parameters. There is no problem with the approximate crossing symme- 
try, even though a solution with a CDD pole can not be computed iteratively 
close to fj,M- 
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We need to include CDD poles in the Pu and P33 partial waves corresponding 
to Delta and Roper resonances. These CDD poles represent an infinite and 
correlated summation of higher order counter terms of the chiral Lagrangian 
Q.For more technical details we refer to |1] . Note that our approach is at vari- 
ance with other approaches where the isobar field is introduced as an explicit 
degree of freedom in the chiral Lagrangian [T5|16|17j . We do not argue against 
an explicit isobar field in the chiral Lagrangian. In fact a possible extension of 
our approach would include an explicit isobar field in the chiral Lagrangian. 
However, so far all our results indicate that the most important ingredients for 
successful descriptions and predictions of pion and photon induced reactions 
off the nucleon based on the chiral Lagrangian are unitarity and analyticity. 



3 Spin polarizabilities of the proton 



The spin structure of the nucleon is described by the effective interaction of 
the third order 



= -47r 



ieff = -47r ^1 -fElEl a-{Ex E) + \ -fMlMl a ■ {H X H) 

- \ 1MIE2 {f-aEi + ViE ■a)Hi + \ ^eim2 {V ■ a + ViH ■ a) E^ 



(5) 



where the four spin polarizabilities 7_biei,7mimi, 7mi_B2, 1eim2 are related 
to the multipole expansion [18]. Spin polarizabilities can be also defined in 
terms of threshold values of invariant scattering amplitudes. Here we apply 
a set of amplitudes, F^_^{y/s,t), introduced in |i4j. They are particularly well 
suited for the construction of partial-wave amplitudes, T^atiV^)^ ^^at are 
free from kinematical constraints and therefore useful in dispersion-integral 
representations like ([2]). Owing to the MacDowel relations 



+^s,t) = Fr'{-Vs,t), (6) 

there are in fact only three independent amplitudes [19|i4j . Similar MacDowell 
relations hold for the partial-wave amplitudes 



7 ^ I +TI „J—^/s) for a = b 

r-^.(+v^)= ^: , (7) 

-Tla,i-V~s) for a^b 

where a,b = 1,2 run over the two possible photon spins. Note that the more 
conventional multipole amplitudes introduced in [TH] do not satisfy the Mac- 
Dowell relations ([T]). The Compton scattering tensor T^,^(g, g; w) reads 
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T,u{q. q;w)=Y.(FtiV:s, t) iJ^f^ {q, q, w) + Fr(^s, t) iJ^'^ (g, g, w) 
i=i ^ 




where q and q are the 4-momenta of the initial and final photons respectively, 
and w is the total 4-momentum of the photon-nucleon system. 

The spin polarizabilities are most economically identified in terms of threshold 
values of amplitudes, Ai_q{v, t), that are even under a crossing transformation 
[20] . We derive their relation with the MacDowel amplitudes 



s — zt 

Ai{-u, t) = +A{+u, t) , v=- , 

1 , , , m,^, — 2 mr,\/s — 3 s — 2 z/ y/s . ^ ^ 
= 2 {Ft + F+) + ^^^-= ^^3+ + (v^ ^ -y^) , 
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TTl 
Z V \J s 



The amplitudes Aili', t) are free of kinematic singularities and constraints and 
are even functions of ly, which is a consequence of crossing symmetry. The spin 
polarizabilities are expressed in terms of the coefficients a^, that are defined 
by threshold values of Born-subtracted amplitudes, y4i(z/, t). It holds 



a, = A,(0,0)-Af°™(0,0), 

a2 — + 2a5 + ag a2 — — ag 

lElEl = t:. , 7-B1M2 = t:. 

Sir rUp OTT nip 



7 



_ 02 + 04 + 2a5 — fle _ ^2 + ^4 + , . 

IMIMI — ^ ' '1M1E2 — n • y^^J 



in nip svr nip 



While our approach respects the unitarity and causahty constraints exactly, 
crossing symmetry is conserved at an approximate level only. As a consequence 
the amplitudes Aj(z/, t) will not be even functions in v. Depending on the 
detailed form of the crossing violations some amplitudes could even be singular 
in the limit z/ — 0. This is not the case in our approach and therefore it is 
possible to identify the spin polarizabilities in an unambiguous manner. 

From ([9]) it follows that the regularity of the amplitude Ai{y^ t) for i = 3, 4, 5, 6 
requests the relations 

F±K,0) = F2±K,0), (11) 

for the Born subtracted amplitudes. This is realized in our approach due to 
the matching condition at a/s = nip in Q that insures that the partial-wave 
amplitudes at this point coincide with their xPT values [1]. 



In the Appendix |b] we express the invariant amplitude Ff[^,t) in terms of 
the partial- wave amplitudes T^^^{y/s). As an application we find 



d 



7 El El 



IQn nii d y/s 



^ATl^,,{nip)-ATln{mp 



+ 



327rm4 



AT 



7M1M1 



IGtt d a/s 



AT 



+ 



4 3 

-^ATl^iirrip) 



AT_2„(mp) 

AT_^,n(mp) 
ATi^nK)' 



7E1M2 



7Afl_B2 



327rm4 



327rm4 



ATVK)-^AT+^ 



,21V^P 



(12) 



where ATj^^^(y^) stays for the non-Born part of the partial-wave amplitude. 
Our result (12) has an interesting consequence. While for the polarizabilities 
'jEiM2,'jMiE2 we recover the results of strict xPT, for 'Jeiei,'Jmimi rescat- 
tering effects are involved in terms of threshold derivatives of partial-wave 
amplitudes. 

At present there are empirical constraints on two combinations of the polariz- 

j. The forward (70) and backward (7,^) polarizabilities 



abilities 
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are 



70 — —lElEl — IMIMI — 1EIM2 — 1M1E2 

= (-1.01 ± 0.08 ± 0.13) 10" W , 



(13) 



Ttt — —lElEl + IMIMI — 1E1M2 + 7M1E2 

= (8.0 ± 1.8) 10~ W . 



(14) 



The comparison of our results for the spin polarizabihties with other theoret- 
ical approaches is presented in Table [l] Our results for the forward and back- 
ward spin polarizabihties are in agreement with experimental values within 
error bars. 

The results from strict xPT in the second and third columns of Table. [T] 
follow from the formal results presented in [3] but using the parameter set 
of |4j. Comparing the with the results is significant since both results 
are determined by the same parameters [3]. As can be seen from Table, [l] 
strict xPT appears not convergent for ""yEiEi and ■Jmimi- In contrast for the 
off diagonal elements ■Jeim2 and 'Jmie2 one may hope for slow convergence. 
This is interesting in view of the large rescattering effects obtained for ■Jeiei 
and jMiMi in our approach. While for 'Jeim2 and 7MIE2 we recover the xPT 
results, for 'Jeiei and 7mimi we predict significant modifications from the 
xPT results. It is amusing to see that our results for •Jeiei and 7mimi are 
quite compatible with the prediction of [26j based on fixed-t unsubtracted 
DRs. However, note some discrepancies for •jEiAn and 'yMiE2- 

The convergence pattern for •Jeiei and 7mimi is reminiscent of the one for 
the s-wave tt" photoproduction multipoles There are indications that in 
our scheme there is significant accelerated convergence as compared to the 
xPT convergence. In order to consolidate this expectation one may perform 
computations within our method to higher orders in the chiral Lagrangian. 





XPT, |3] 


xPT, m 


DK |26|27| 


this work 


lElEl 


-5.93 


-1.41 


-4.3 


-3.68 




-1.19 


3.38 


2.9 


2.47 


1E1M2 


1.19 


0.23 


0.0 


1.19 


1MIE2 


1.19 


1.82 


2.1 


1.19 


70 


4.74 


-4.02 


-0.7 


-1.16 


1-K 


4.74 


6.39 


9.3 


6.14 



Table 1 

Proton spin 
polarizabilities 
obtained from 
different sources 
in units of 10""^ 
fm^. 
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Fig. 1. The beam asymmetry S3 for linearly polarized photon and unpolarized 
proton target, as function of photon laboratory energy and for different center of 
mass scattering angles as indicated. The solid curves result from our chiral approach. 
The dashed and dotted curves follow from subtracted DR calculation as explained 
in the text. The data are from [28j and correspond to the three scattering angles 
shown in the figure. 

4 Spin polarization in Compton scattering 



The beam asymmetry S3, beam-target asymmetries Ti2z and Ti2x for the proton 
Compton scattering will be measured at MAMI [7j in a wide energy range, 
starting from threshold. We refer to the Appendix |X] for their definitions. 

The main purpose of these measurements is the determination of the spin 
polarizabilities of the proton. The appropriate tool to analyze the data set 
are subtracted fixed-t DRs [25]ITU] . which use as input empirical pion photon- 
production amplitudes. In this approach the spin polarizabilities are free pa- 
rameters that will be adjusted to the data set. In order to shed further light 
on the intrinsic theoretical uncertainties it may be useful to compare with our 
present approach based on the chiral Lagrangian. 

We recall the definitions of the various spin observables in terms of helicity 
matrix elements, ^^(a/s,^), and helicity partial-wave amplitudes, tj_^fj{y/s), 
in the convention used in The helicity matrix elements are expressed in 
terms of the helicity partial-wave amplitudes as follows 



1 COS I 



P' 1 (cos 

2 
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AtT ^/s 




-^)(J+i) 










sin 6 cos 1 


471 y/s 


V ^ 








-Pj- 


Trip 


2 


cos^ 1 ^ 






l)(J+i) ^ 


+ iJ 


-1) 




nip 




2 sin^ 1 




i (J 


-i)(J+|) 


+ {J 




p;^i(cos0) 



}. 



1 (^cos 

2 



E«22-<22) {3(J-^)P;_,(cOSe) 

+ 2 (P;: 1 (cos ^) + P;: 3 (cos e))], (15) 



with the scattering angle 9 in the center of mass frame. The spin observables 
have a transparent representation in terms of the hehcity matrix elements 



E3 ^ = 3ft ((01 + 05) 03 + (02 - 06) 0:) 



+ 



+ 



-'2x 



da 1 



52 + 06)<^3j , 
^l|'-|0 



(16) 



'dn 2 

For the readers convenience we recall the relation of the helicity partial-wave 
amplitudes with the covariant partial- wave amplitudes used in (12). It holds 




Pcra 

2 J-l mp 

2J+3 Pen 







(17) 



db 



Since in our approach there are no free parameters we compare our results 
(solid curves) in Figs. [T]j3] to calculations based on subtracted fixed-t DRs 
(dashed and dotted curves) . For the dashed curves we use as input our predic- 
tions for both the spin-dependent and the spin-independent polarizabilities. 
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Fig. 2. Double polarization asymmetry Y^2x for circularly polarized photon and tar- 
get spin aligned transversely w.r.t. the photon momentum, as function of laboratory 
photon energy and scattering angles as indicated. The solid curves result from our 
chiral approach. The dashed and dotted curves follow from subtracted DR calcula- 
tion as explained in the text. 




[GeV] 

Fig. 3. Double polarization asymmetry Ti2z for circularly polarized photon and tar- 
get spin aligned longitudinally w.r.t. the photon momentum, as function of labo- 
ratory photon energy and scattering angles as indicated. The solid curves result 
from our chiral approach. The dashed and dotted curves follow from subtracted DR 
calculation as explained in the text. 

The dotted curves are obtained with the empirical spin-independent polariz- 
abilities [22] together with the values of the spin polarizabilities given in the 
second last column of Tab. 1. The difference between our predictions and the 
DR calculation can provide a useful cross check and estimate residual uncer- 
tainties. For a discussion about possible sources of uncertainties in the DR 
calculation we refer to ITOl. 



In Fig. [T]the results for the beam asymmetry are compared with available data 
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j28] at different scattering angles. New measurements at MAMI are expected 
to have a much better accuracy, typically at a few percent level. In Figs. |2]j3]we 
show the results for the double polarization observables Y^^x and Y^^z which are 
mostly sensitive to the 7eiei and 7mimi polarizabilities, respectively. The solid 
and dashed curves agree well below the pion production threshold and start 
deviating slightly when the energy increases. The small difference between the 
two curves illustrates the precision one may expect for our results. Though the 
two methods are quite different in construction, they both rely on analyticity 
and unitarity which impose stringent constraints in the results for Compton 
observables. The closeness of the dashed and dotted curves emphasize the 
similarity of the polarizabilities as predicted by our chiral and unsubtracted 
DR approach. 



5 Summary 

We studied nucleon Compton scattering from threshold up to Eiah = 350 MeV 
within an approach developed in |3] and based on an analytic extrapolation of 
subthreshold amplitudes calculated in xPT. The Compton sector involves no 
additional free parameters. Predictions for the spin polarizabilities were made 
within this scheme. We found a nice agreement of the predicted values for 
the forward and backward spin polarizabilities 70 and 77r with experimental 
values. All four spin polarizabilities 'Jeiei, 7mimi, 1eim2, 7mie2 are in a rea- 
sonable agreement with the calculations based on fixed-t DRs [26p0] . While 
for the transition polarizabilities 'yEiM2,lMiE2 "we recover the results of strict 
xPT, for the diagonal 'Jeiei, Imimi elements we find significant effects from 
rescattering. 

We also analyzed our predictions for several spin observables that will be 
measured at MAMI [7J. Our results are close to the ones obtained from fixed- 
t DRs, although the two approaches are quite different in construction. This 
illustrates again the enormous predictive power of combining the unitarity and 
analyticity constraint in hadron physics. Furthermore, the comparison between 
the two calculations provides an important cross check and can be used to 
estimate the theoretical uncertainty in calculations of Compton scattering 
observables. 
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A Appendix 



The beam asymmetry for photons which are linearly polarized either parallel 
or perpendicular to the scattering plane and unpolarized nucleons is defined 

as 



where a denotes the yield of particles corresponding to a particular polariza- 
tion of the initial photon. 

The asymmetries with circular photon polarizations (left-hand or right-hand) 
are defined as 



= ' = (A.2) 

for the nucleon spin lying in the reaction plane and aligned in ±x or ±z 
directions, respectively. For more details we refer the reader to Ref. |T8] . 



B Appendix 

We decompose Ff_^[^^,t) in terms of partial wave amplitudes 
where we detail all non vanishing coefficients with 



J,m,n '"'P 



J 



rrip 2cjw^ 2 {2J + 3)w^ " ''^^ 

nip 8cjw^ +'2^ 8(2J + 3)w7 +-22 

(l-u;f(l + t^)2 (1 + 2^-^2) ^^ 
+ 8cjw^ 
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+ 



{l-3w){l-wf{l + w) {l + 2w-w^) 

8 (2 J + 3) w7 



,22 ) 



rrir, 



2\5 



h,22 



,22 ) 



Ff (1 



2\7 



32C>15 ^+,22 



+ 



(l-w)^(l + w)5(l + 2w-w2) 



32 c5 wis 



-,22 ' 



(B.2) 



with cj = a/(2 J — 1) (2 J + 3) and w = ^fslnip and 



_ / 2J-1 1-.^ , 1-2J , ^T^J 

m, "V2J + 3 2«;2 ^+'2i+2(2J + 3)«; +'^^^4«;3 -'11 



/2J-1 l + w^ 1-2J 
2J + 3 2^2 -'^^ ~ 2(2J + 3)«; 



'^ = ^^^(6J + l + (2J + 3) (u;^ + ^;-1)w)t^2i- "^7 



5 - 2 J + (2 J + 3) (3 tt;-'^ + 4 u;^ - 2 w - 4) . 

8(2J + 3)w5 +'22 
/ 2J + 3 (1-^)3(1 + ^)2 (i_3^)(i_^)2(i + ^) 

V2J-1 8«;6 -'21+ 8^5 ^-,22, 



F20_ -1-2^3 ^3^4 1-^3 J 

2cji/;6 ^+'21+ (2J + 3)«;4'+'22 
+ 2cj^/;6 -'21+ (2J + 3)^/;4 ' 

F21 _ (1 - ^^;2)3 (1 + ^2) _ (1-^2)2 



m. 



(^2J + 5-8(J+l)w 



,22 



+4 (1 - 2 J) ti;2 + (2 J - 7) (w - 2) j T 

(1 - wf{\ + w)2 (1 + w + w2 + 3^2) J 
+ 8^7^;;^^^ -'21 

_ (l-^)2(l + w) (2J + 5-3(2J+l)w-4(4J+l)w2 

+2 (7 - 2 J) + 9 (2 J + 1) + (7 - 2 J) v/^ , 

Ff _ (1 - ^2)4 

mp ~ 32 (2 J - 1) u;i3 
(l-t^)7(i+^)5 

^ 32 (2J- 1)^13 -'22' 

Ff (1-W)5(1 + W)3_; (1-^2)2. „ , 

— = ^ 7^2^ -^+22+ \ 2 Q (1-2W + 2W' + W^)T:^22, 
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J 

22 
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;i + / , o 2 , . 3 4 , 5\nnJ 

1 + w + 2w + 4:W — w + w ]T_ 



{l-w)\l+wY 

(l-wfil + w)' 

{l-w)%l+wY 
64 4 



22 



1 - w + 3 - T 



l + w + 3w^-w^ T_ 



,22 ) 



,22 



l + 3w + 6w' + 6W^ -5w^ + w^]T 



,22 ) 



(B.3) 



and 
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2(1-^) w 2 

' +'2^ (2J + 3)«;3 
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T 



J 



T 
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22 -^-,22 
2\ (Q A „„ Q„„2\ 



Cj W 



3 ^-,21 ) 



(l-t^^) (3-4u;-3^^ ^^ 

2cjw' +'21 2(2J + 3)u;7 '+'^2 

ujfil+wY {l-3w){l-wY{l+w) 
2cjw^ -'21 2(2J + 3)^7 ^-,22, 



{l-wf{l+wf J {l-wy{l-2w-w^) 



rno _|_ 

' ,22 
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2 4u'ii 



15 



2 4^11 
(1 -uO'fl 



rpJ 

^ -,22 ' 



.22 



CjW 
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(B.4) 



To obtain the expression for one can use the identities F^{—w,t) 
Fr{w,t) and (B. Note that T^^l^'^ = for a = 2 or 6 = 2. 
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